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Abstract. Modified f{R) gravity in the Palatini approach has been presently applied 
to Cosmology as a realistic alternative to dark energy. In this concern, a number of 
authors have searched for observational constraints on several f{R) gravity functional 
forms using mainly data of type la supenovae (SNe la). Cosmic Microwave Background 
(CMB) radiation and Large Scale Structure (LSS). In this paper, by considering 
a homogeneous and isotropic flat universe, we use determinations of the Hubble 
function H{z), which are based on differential age method, to place bounds on the free 
parameters of the f{R) = R — /3/i?" functional form. We also combine the H{z) data 
with constraints from Baryon Acoustic Oscillations (BAO) and CMB measurements, 
obtaining ranges of values for n and /3 in agreement with other independent analyses. 
We find that, for some intervals of n and (3, models based on f{R) = R— P/R^ gravity 
in the Palatini approach, unlike the metric formalism, can produce the sequence of 
radiation-dominated, matter-dominated, and accelerating periods without need of dark 
energy. 
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1. Introduction 

Nowadays, one of the key problems at the interface between fundamental physics and 
cosmology is to understand the physical mechanism behind the late-time acceleration 
of the Universe. In principle, this phenomenon may be the result of unknown physical 
processes involving either modifications of gravitation theory or the existence of new 
fields in high energy physics. Although the latter route is most commonly used, which 
gives rise to the idea of a dark energy component (see, e.g., [U El El HI El E] ) , following 
the former, at least two other attractive approaches to this problem can be explored. 
The first one is related to the possible existence of extra dimensions, an idea that links 
cosmic acceleration with the hierarchy problem in high energy physics, and gives rise 
to the so-called brane- world cosmology [3 [HI (U [TOl [11]. The second one, known as 
f{R) gravity, examine the possibility of modifying Einstein's general relativity (GR) 
by adding terms proportional to powers of the Ricci scalar R to the Einstein-Hilbert 
Lagrangian [III [131 [H [I5] . 

The cosmological interest in f{R) gravity dates back at least to the early 1980s and 
arose initially from the fact that these theories may exhibit an early phase of accelerating 
expansion without introducing new degrees of freedom [16]. Recently, f{R) gravity 
began to be thought of as an alternative to dark energy [TTl [181 [El [20] and a number 
of authors have explored their theeoretical and observational consequences also in this 
latter context. As a consequence, many questions have been raised and there is nowadays 
a debate about the viability of such theories (see, e.g., [211 [22l [231 [211 [251 [26] ) . However, 
it seems that most of problems pointed out cannot be generalized for all functional 
forms of f{R)- For example, it has been shown that specific forms of the function f{R) 
may be consistent with both cosmological and solar system-tests [271 [211 [221 [30] . Non- 
linear coupling of matter with f{R) gravity has been discussed by [31], and [321 [33] 
discussed connections with MOND theory as well as comparison with solar observables. 
Besides, by starting from general principles such as the so-called energy conditions, and 
by generalizing them to f{R) gravity, p4] have shown how to place broad constraints 
to any class of f{R) theory. 

Another important aspect worth emphasizing concerns the two different variational 
approaches that may be followed when one works with modified gravity theories, namely, 
the metric and the Palatini formalisms (see, e.g., [35]). In the metric formalism the 
connections are assumed to be the Christoffel symbols and variation of the action is 
taken with respect to the metric, whereas in the Palatini variational approach the metric 
and the affine connections are treated as independent fields and the variation is taken 
with respect to both. In fact, these approaches are equivalents only in the context of 
GR, that is, in the case of linear Hilbert action; for a general f{R) term in the action 
they give different equations of motion. 

For the metric approach, a great difficulty in practice is that the resulting field 
equations are fourth order coupled differential equations which presents quite unpleasant 
behavior. In addition, simplest f{R) gravity models of the type f{R) = R — (3/R^ have 
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shown difficulties in issues such as passing the solar system tests EZ], in having 
the correct Newtonian limit |3H1 139] and gravitational stability [lOj. In a recent study 
[211 [22] have shown that these theories cannot produce a standard matter-dominated 
era followed by an accelerated expansion. 

On the other hand, the Palatini variational approach provides second order 
differential field equations which can also account for the present cosmic acceleration 
without need of dark energy. Recent studies [HI |12] have shown that the above 
cited power-law functional forms are capable of producing the last three phases of 
the cosmological evolution, i.e., radiation-dominated, matter-dominated, and late time 
accelerating phases. Some issues still of debate in literature are whether f{R) theories in 
Palatini formalism satisfy the solar system tests and give the Newtonian approximation 
[441 |45| l46l l47l l48] . and whether they are free of gravitational instabilities [49l l50] . 
Recent reviews on f{R) and other modified gravity theories can be found in [511 [52| [53]. 

From the observational viewpoint, however, it is important to look into whether 
these theories of gravity are indeed compatible with different kinds of currently available 
cosmological data. In particular, the observational viability of some functional forms of 
f{R) gravity have been studied using mainly SNe la, CMB radiation and Large Scale 
Structure (LSS) data [HI [Ml Il2l US] . 

In this paper, by following [55] , we use determinations of the Hubble parameter as a 
function of redshift [56] to derive constraints on the parameters of the f{R) = R — (3/ R^ 
theory of gravity in the Palatini approach. These determinations, based on differential 
age method, relates the Hubble parameter H{z) directly to measurable quantity dt/dz 
and can be achieved from the recently released sample of old passive galaxies from 
Gemini Deep Deep Survey (GDDS) [57], [58] and archival data [59l [60l [61]. The same 
data, along with other age estimates of high-^r objects, were recently used to reconstruct 
the shape and redshift evolution of the dark energy potential ^2], to place bounds 
on holography- inspired dark energy scenarios [63], as well as to impose constraints on 
the dark energy equation of state parameter {w) by transforming the selected GDDS 
observations into lookback time determinations [6l]- We also combine H{z) data with 
BAO [65j and the CMB shift parameters ^66] to better constrain the free parameters of 
our f{R) model. A brief discussion on the cosmic eras in the context of the Palatini 
approach is also included. 

2. Basic equations in the Palatini approach 

The simplest action that defines an f{R) gravity is given by 

S=Ynj d'xV^fiR) + Srn, (1) 

where k = SttG, G is the gravitational constant and Sm is the standard action for the 
matter fields. Here R = g'^^RapiJ'^) and F^^ is the affine connection, which in the 
Palatini approch is different from the Levi-Civita connection 
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By varying the action with respect to the metric components we obtain the field 
equations 

/ 

fRR^,u{T) - -g^,u = nT^u , (2) 

where fn = df/dR and T^^, is the matter energy- momentum tensor which, for a perfect- 
fiuid, is given by T^^, = {pm + Pm)u^Uy + Pmg^iv, where pm is the energy density, pm is 
the fiuid pressure and is the fluid four-velocity. Variation of action ([T]) with respect 
to the connection provides the equation that determines the generalized connection: 
V/3 [/i?V~fi'fi''"^] — O5 where V is the covariant derivative with respect to the affine 
connection F^^. This equation implies that one can write the conection F as the Levi- 
Civita connection of a conformal metric 7^1, = fng^v [671 EH]- The generalized Ricci 
tensor is written in terms of this connection as 

p (J^\ "p« -pa _|_ -pa -pA -po! pA /o\ 

We next consider an homogeneous and isotropic universe and investigate the 
cosmological dynamics of f{R) gravity in a flat Friedmann-Lemaitre-Roberston- Walker 
(FLRW) background metric g^y = diag{—l,a'^,a'^,a'^), where a{t) is the cosmological 
scale factor. By expressing the generalized Ricci tensor ([3]) in terms of the Ricci tensor 
R^y{g) associated with the metric g^i, we obtain the generalized Friedmann equation 
[69] 




/ = I^Pm (4) 

where H = d/a is the Hubble parameter and a dot denotes derivative with respect to 
the cosmic time t. Here, we adopt the notation fji = df/dR, fjm = d'^f/dR^ and so on. 
The trace of Eq. ([2]) gives 

fRR-2f = -KPm, (5) 

where we have considered the fluid as a pressureless dust satisfying the conservation 
equation p^ + ^Hpm = 0. By combining this equation with the time derivative of Eq. 
([5]) we flnd 

R = ^^^. (6) 

JRRJ^ — JR 

Now, by substituing Eq. into Eq. (jlj) we obtain 

+ f^ R - f 

QfRi 
where 

^ = 1 _ 3 fRRjfRR - 2/) 
2 fRifRRR — fR) 

Note that the usual Friedmann equations are fully recovered from the above expressions 
if f{R) = R, in which case the action ([1]) reduces to the Einstein-Hilbert one. 



2 ^l^Pm^ JRJ^- J /„N 

^ = FTT^ ' V) 
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2.1. Parameterization 

In this work we are particularly interested in testing the viability of a general functional 
form given by 

fiR) = R-P/R-. (9) 

In a recent paper, [H] have found that this model can be compatible with the supernova 
"Gold" data set from [70] for a given interval of the parameters (3 and n. More recently, 
[12] have shown that models of this kind are compatible with the Supernova Legacy 
Survey (SNLS) data [71] and also found narrow ranges for the values of n and (3 
consistent with that from [51]. Here we will follow the numerical scheme used by [12] 
to obtain the Hubble function H{z). 

Firstly, we rewrite Eqs.Q and ([7]) in terms of redshift parameter z = ao/a — 1 and 
the density pm = Pmo(l + zf- 

fnR - 2/ = -?,Hln^o{l + zf , (10) 

and 

_ 'inmo{l+zf + f/Hl 

Hi 6fRe ^ ' 

with 

2 fn RfuB^ - Ir 

where VLmo = i^pmo/ {"^Hl). An important aspect worth emphasizing at this point is 
that Eqs. (fTOl) and ( fTTI) evaluated at z = impose a relation among n, VLmo and [3, so 
that specifying the values of two of these parameters the third is automatically fixed. 
In other words, in the Palatini approach, a f{R) = R — j3/R^ theory introduces only 
one new parameter: n or (3. In the following, we will always work with n as the free 
parameter. 

3. Analyses and discussion 

In order to impose constraints on models of f{R) gravity given by Eq. (Q, we minimize 
the function 

2 _ [HthjZils) - Hobs{,Zi)f 

k ^'(^^) 

where Hth{zi\s) is the theoretical Hubble parameter at redshift zi given by ffTTl) which 
depends on the complete set of parameters s = {Ho^VLmo^n); Hobs{zi) are the values 
of the Hubble parameter obtained from the data selected by [62] (SVJ05) and <j{zi) is 
the uncertainty for each of the nine determinations of H{z). In what follows, we will 
work with n and flmo as free parameters and study the bounds on them imposed by 
the SVJ05 H{z) data sample. We also marginalize over the current value of the Hubble 
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Figure 1. Left: The Hubble parameter H{z) as a function of the redshift for the best- 
fit values of n and 0,mo using H(z) data only and a combined fit including BAO and 
CMB shift measurements. For the sake of comparison, the standard ACDM model 
prediction is also shown. The data points are the measurements of the H{z) by 
[52]. Middle: Contour plots at 95.4% and 99.73% c.l. in the n x n„,o plane for 
a f{R) = i? — (3/ R"" theory using the SVJ05 sample of H{z) measurements. Right: 
Effective equation of state [Eq. p^ ] as a function of redshift for the best-fit value of 
n and f2mo from H{z) data analysis. A radiation component with f2^o = 5 x 10~^ has 
been included. 



parameter by assuming the Gaussian prior Hq = 72 ± 8 km/s/Mpc, in agreement with 
the final results of the HST key project [72]. 

In Figure la we show the evolution of the Hubble parameter with redshift for the 
two best-fit values for n and fi^o discussed in this paper, as well as the prediction from 
the standard ACDM model (fimo = 0.27). The three curves are superimposed on the 
data points of the SVJ05 sample. Note that all models seem to be able to reproduce 
fairly well the H{z) measurements. 

Figure lb shows the first results of our statistical analyses. Contour plots (95.4% 
and 99.7% c.l.) in the n x VL^^ plane are shown for the given by Eq. f|T3l) . We clearly 
see that the measurements of H{z) alone do not tightly constrain the values of n and 
VLmoi allowing for a large interval of values for these parameters, with n ranging from 
-1 to even beyond 1, and VL^o consistent with both vacuum solutions (fimo = 0), as 
well with universes with up to 70% of its energy density in the form of nonrelativistic 
matter. The best-fit values for this analysis are = 0.04, n = —0.90 and (3 = 1.11. 

Figure Ic shows the effective equation of state 
^ 2{1 + z)dH 

-^^^ = -^ + T^ 

as a function of the redshift for the best-fit values above. To plot this curve we have 



included a component of radiation Vt^o = 5 x 10~ . It is worth mentioning that the 
best-fit point is not representative from the statistical point of view, given the weak 
power of constraining shown in the Figure. Note also that, similarly to some results 
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Figure 2. Left: Same as in Fig. 2 when BAO and CMB shift parameters are 
included in the analysis. Middle: Same as in Panel 3a when BAO and CMB shift 
parameters are included in the statistical analysis. Right: Effective equation of state 
as a function of redshift for a f{R) = R — /3/i?" theory with n = —0.15 and different 
values of flmo — 0, 0.01, 0.05, 0.25. Note that a direct transition from radiation to a de 
Sitter era happens as flmo approaches zero. 



in the metric formalism [211 122] , for these specific values of the n and 
there is no matter-dominated era followed by an accelerated expansion. 



parameters, 



3.1. Joint Analysis 



In |l2] it was shown that when the measurements of SNe la luminosity distances are 
combined with information concerning the Baryon Acoustic Oscillation (BAO) peak 
(measured from the correlation function of luminous red galaxies) and the CMB shift 
parameter (which relates the angular diameter distance to the last scattering surface 
with the angular scale of the first acoustic peak in the CMB power spectrum), the 
constraining power of the fit to f{R) paramenters is greatly improved. Following such 
an approach we examine here the effects of summing up the contributions of these two 
parameters into the of Eq. f|T3l) . 



In fact, when the BAO |65| 



A 



2bao=0.35 

^BAO 

and the CMB shift parameter | 



/72 





/ ^BAO \ 
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1 




^BAO 




\i ") 











1/3 



0.469 ±0.017, (15) 



1089 



R 



1089 



dz 



1.70 ±0.03 



(16) 



X To include the CMB shift parameter into the analysis, the equations of motion had to be integrated 
up to the matter/radiation decoupling (z ~ 1089), so that radiation is no longer negligible and was 
properly taken into account. 
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this paper 
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this paper 


0.03 


4.7 



Table 1. Comparison with independent determinations of best- fit values for n and (5 
(the ACDM model corresponds to n = and /? = 4.38). 

are included into the fit, a considerable enhancement of the constraining power over n 
and f2mo takes place, as can be seen in Fig. (2a), which shows the contour curves in the 
n X fi„o plane. The best-fit value (n = 0.03, /? = 4.7, fi^o = 0.26 with xV^dof ^ 1.02) 
is consistent with current estimates of the contribution of non-relativistic matter to the 
total energy density in a flat universe. The fit also constrains the parameters n to lie in 
the intervals n G [—0.25,0.35] and /3 G [2.3,7.1] at 99.7% c.l., which is consistent with 
the results obtained in Refs. [HI using the supernova Gold and the SNLS data sets, 
respectively. For the best-fit solution, the universe goes through the last three phases of 
cosmological evolution, i.e., radiation-dominated [w = 1/3), matter-dominated {w = 0) 
and the late time acceleration phase (in this case with w ~ —1), as shown in Fig. 
(2b). For the sake of completeness and also to better understand the role of Qmo on the 
behavior of the effective equation of state in an /(i?)-theory, we plot in Fig. (2c) Weff 
as a function of z for 4 different values of ^mo, namely, 0.0, 0.01, 0.05, 0.25 for a fixed 
value of n = 0.03. As physically expected, in the limit flmo = we clearly see that the 
matter era at intermediate redshifts is completely absent. 

3.2. Future data 

In [62], the constraining power on the dark energy potential from the future Atacama 
Cosmology Telescope (ACT) data [73] was estimated. Such a project is expected to 
identify ~ 2000 passively evolving galaxies, hosted by clusters with masses greater than 
IO^^Mq. The clusters will be identified through their Sunyaev-Zeldovich imprint in the 
CMB, and spectroscopic information of their brightest member galaxies obtained by 
telescopes in South Africa and Chile. Assuming that ACT will be able to obtain ~ 
1000 measurements of the Hubble parameter, the authors were able to make forecasts 
in the phase space of dark energy potential parameters. Here, we will be a little more 
conservative and will assume that ~ 200 determinations of H{z) will be available in the 
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m 

Figure 3. Predicted constraints in the n x 0„io plane (99.7% C.L.), for a sample of 
200 determinations of H{z) expected from the Atacama Cosmology Telescope (ACT) 
[73] in the redshift range 0.1 < z < 1.5. A 15% error in H{z) has been assumed. The 
bulk model is taken to be consistent with the ACDM model for ^Imo = 0.27. 

interval 0.1 < z < 1.5. As in [62], we will assume a 15% uncertainty in each value. In 
order to simulate the future n x flmo parameter space, we have adopted a bulk model 
consistent with the ACDM model with fi^o = 0.27 for the H{z) low redshift values. 

Figure [3] shows the final constraints for such a future data set with and without 
combining it with the BAO and CMB shift parameter. We clearly see that low redshift 
H{z) determinations alone do not constrain tightly the /(-R)-theory parameters. Note 
also that, when combined with high- 2; information, such as the one provided by the 
CMB shift parameter, the constraining power is highly enhanced for the combined data 
set, but the difference compared to the current fit is not that big. In Table I we 
summarize the main results of this paper compare them with recent determinations of 
the parameters n and f3 from independent analyses. 

4. Conclusions 

By considering a flat FRW cosmology we have analyzed the f{R) = R — theory of 
gravity, with equations of motion derived from the generalized action ([T]), according 
to the Palatini approach. We have performed consistency checks and tested the 
cosmological viability for a theory of this type by using current determinations of the 
Hubble parameter at different redshifts obtained from differential ages techniques. The 
use of these data to constrain cosmological models is interesting because, differently 
from distance measurements, the Hubble parameter is not integrated over. This means 
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that the differential age method is less sensitive to systematic errors than the standard 
distance methods. We find that the determinations of H{z), when combined with the 
BAO and CMB shift parameter, lead to constraints competitive to those achieved with 
SNe la Gold and SNLS data, as given by [HlllS]. The FRW cosmology corresponding 
to the best-fit solution for a combined iJ(z)+BAO+CMB minimization presents all 
three last phases of the Universe evolution, namely, radiation era, matter era and a 
phase of acceleration at late times. We emphasize that a great improvement on the 
H{z) measurements is necessary for the differential age technique to provide strong 
additional constraints on cosmological parameters such as those coming from SNe la, 
CMB and LSS data. As we have shown, the combination of 200 determinations of 
H[z) at a 15% accuracy level (as expected from the ACT project) with BAO and CMB 
shift parameters was shown to provide constraints on f{R) variables similar to the ones 
already obtained with the current H{z) measurements. 
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